gusucode.com > 《MATLAB图像与视频处理实用案例详解》代码 > 《MATLAB图像与视频处理实用案例详解》代码/第 14 章 基于主成分分析的图像压缩和重建/pcasample.m
function [coeff,score,rate]=pcasample(X,p) [V,D]=eig(X'*X); for i=1:size(V,2) [~,idx]=max(abs(V(:,i))); V(:,i)=V(:,i)*sign(V(idx,i)); end [lambda,locs]=sort(diag(D),'descend'); V=V(:,locs); coeff=V(:,1:p); score=X*V(:,1:p); rate=sum(lambda(1:p))/sum(lambda);